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The subject of this report is the dynamics of elastic system in conjunction with hydrogen bonds 
of the DNA. We draw attention to the draw-back of the familiar rod model of the DNA, and make 
a case of constructing models that could accommodate the intrinsic structure of the DNA. In this 
respect studying the interplay among the elastic system and the protons of the DNA, is of interest, 
for it could accommodate the inter-strand as well as the tunneling modes of protons. Following 
this direction, we come to the conclusion that the elastic-proton dynamics may have a bearing on 
biophysics of the DNA. The phenomenon of point mutations is discussed within this framework. 

PACS numbers: 87.15- v 

Keywords: DNA, hydrogen bonds, proton tunneling, mutations 



INTRODUCTION 



A molecule of the DNA can attain several hundred /im in length. If we neglect details that have a 
size of one thousand A, or less, we can visualize it as a soft shapeless line and conclude that on this 
scale it behaves like an ordinary polymer. In contrast, looking at its smaller segments, of one hundred 
A or less, we observe that it tends to be straight. Thus, borrowing a comparison from everyday life, we 
may say that a molecule of the DNA looks like a piece of steel wire whose long segments are flexible and 
the short ones are stiff. The elastic properties of the DNA are intimately related to its being a double 
helix. The latter imposes severe constraints on deformations which can be effected without destroying the 
molecule and to a large extent determines its mechanical properties. In fact, the two strands comprising 
the molecule of DNA have just small bending rigidities, just as usual polymers. But the formation of 
the two-stranded structure drastically changes the DNA by making it both stiff and capable of forming 
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sophisticated spatial shapes. Similarly, the dynamics of the DNA is largely determined by the relative 
motion of its strands. The functioning of the DNA also involves the charge transport in the molecule. The 
latter may be related either to the motion of electrons along the strands, or to the tunneling of protons 
inside the hydrogen bonds connecting the bases. The problem that still waits its answer is whether there 
is any interaction between the elastic dynamics of the double helix and the tunneling dynamics of charge. 

The diameter of the DNA is about 20A, and taking into account its enormous length we may consider 
it both a microscopical and macroscopical object. Thus, the DNA requires a special means for analysis of 
its physical properties. Important, the approach from 'first principles' aimed at conformational dynamics 
, for example using the methods of quantum chemistry, often comes across very serious difficulties, 
which do not have just technical character, but correspond to the need for a physical picture to rely 
upon. Therefore, the theoretical study of the DNA is to utilize models, which are necessarily based 
on extremely crude simplifications. Their connection with the microscopical structure of the DNA, is 
a matter of intuition rather than of rigorous demonstration. Yet it would not be quite appropriate to 
accept only rigorously demonstrated facts and to ignore ideas that could motivate the development of 
the DNA theory and provide a stimulus to further experimentation. 

Owing to the large stiffness of the molecule of DNA, described by the persistence length of about 50QA 
at which it strongly resists strains caused by heat fluctuations, one may try to visualize it as an elastic 
rod. The above model serves the basis for a number of theoretical approaches, but even though it has 
turned out to be successful, for example for the topological analysis of the DNA conformations, there 
are situations in which it does not work properly; specifically, when external charges are to be taken into 
account. It is also important that it does not allow for the internal degrees of freedom due to the two- 
stranded helical structure. The latter imposes stringent constraints on possible strains of a molecule of 
the DNA, which result in the specific structure of its vibrational modes, very important for understanding 
its functioning. Therefore, we shall consider consequences that may be inferred from the basic properties 
of the DNA, aiming at a qualitative approach that uses simple theoretical models that may accommodate 
the internal degrees of freedom. At this point we should like to note that since different problems require 
the use of appropriate approximations and specific choice of dynamical variables, there is no unique form 
of the energy functional, even though the object of our study remains a molecular of the DNA. Indeed, 
we shall choose various forms of the energy owing to the necessity to accommodate conformations of the 
DNA under consideration and external conditions in action. 

Let us recall that the double helix of DNA consists of long chains, or strands, which have the backbones 
composed of sugar and phosphate residues, and special chemicals, bases, keeping the two strands together. 
The fundamental building blocks of the strands are nucleotides, joined to each other in polynucleotide 
chains. The nucleotide consists of a phosphate joined to a sugar (2'-deoxyribose), to which a base is 
attached. The sugar and base alone are called a nucleoside. The chains, or strands, of the DNA wind 
round each other in a spiral forming a double helix, the bases being arranged in pairs: adenine - thymine 
(AT), guanine - cytosine (GC), so that the sequence of bases in one strand determines the complimentary 
sequence of bases in the other and constitutes the genetic code stored by the molecule of DNA. There are 
several forms of the DNA, denoted by A,B, and Z. The most common one in nature, is the so-called B- 
form. One turn of the helix of the B-form, corresponds approximately to 10.5 base-pairs, and the distance 
between adjacent pairs of bases is approximately 3.4A In real life there are considerable deviations from 
the canonical B-form of the DNA. Therefore, there is a need for a special nomenclature for describing 
its conformations (see [l| for the details). In fact, it is easy to see that generally a considerable set of 
parameters are required to this end. For example, even if we assume a simple picture of base pairs as flat 
plates, we still need to use vector quantities for describing their distance from each other and angles for 
mutual orientation, that is on the whole nine parameters. The most common one, Roll and Tilt, shall 
describe the planes of successive base pairs not being parallel, as is prescribed for the B-form. In fact, 
it is necessary to introduce the so-called propeller angle for defining the measure for deviations of the 
normals to the bases inside a base-pair. It is worth noting that the deviations from the canonical form 
are by no means small, and may have a size of tens of degrees. 

In contrast to the sugar-phosphate backbone of a single strand of DNA, which is formed by strong 
covalent forces, the double helix of the DNA is due to the interplay of weak chemical forces. In fact, the 
two chains of the DNA are held together by the hydrogen bonds between complementary bases and the 
stacking interactions between adjacent bases attached to the same sugar-phosphate strand above each 
other in neighbouring pairs, so that on the whole it requires the energy of order lOKcal/mol or several 
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ksT, whereas the covalent bonds are by orders of magnitude larger, see [2j. The main contribution 
to the thermodynamic stability of the double helix is due to the hydrogen bonds which require the 
specific choice of the pairs, adenine - thymine, guanine-cytosine. The important point about it is that 
there is a competition between the hydrogen bonds formed by base-pairs and the hydrogen bonds of 
the molecule with water in aqueous solution. It is generally believed, even though it is worthwhile to 
mention that there are no adequate theoretical estimates to the effect (see p fTBI where we discuss a rough 
qualitative model of the phenomenon) that the competition makes for the increase of entropy, and thus 
the stability of the DNA. The stacking interaction is due to the bases being flat water-insoluble molecules, 
lying above each other roughly perpendicular to the direction of the helical axis so as to enable electron 
clouds between bases to contribute to the helical stability. This is again only a plausible hypothesis, or 
even simply an intuitive picture that needs careful investigating. The energy of the interaction between 
complementary base pairs has been estimated by various means which rely on experimental data and 
computer calculations within the framework of the quantum chemistry, Q . It is generally falls within 
lOKcal/mol (see Q for more recent results). 

The large persistence length of the DNA, ss 500A, which is more than 20 times larger its diameter, 
serves the main, and essentially intuitive, argument for considering it as an elastic rod and employing the 
methods of continuum mechanics for its study. Strictly speaking, this is not correct, for the cross section 
of this 'rod' corresponding just to a base-pair comprises only several tens of atoms. The main argument 
in favour of such approach is generally the pragmatic one, "the truth is useful" . In fact, as follows from 
numerical simulation within the framework of this approach, the above model appears to be acceptable for 
many regimes that involve functioning of the DNA , [H , @ • It needs some modifications so as to take into 
account the important effects of the anisotropy due to unsymmetrical positions of constituent nucleoside, 
and thus requires the use of anisotropic elastic modulii. Equally important, there is also a geometrical 
asymmetry generated by the relative positions of nucleosides inside a base pair. The asymmetry can be 
visualized as two groves on the surface of the rod. The grooves are helpful in describing the interactions 
of the DNA with external charges by allowing the graphic representation of their distribution on the 
surface of a molecule of DNA. It should be noted as well that a molecule of the DNA is itself negatively 
charged. Therefore, the emerging picture of the DNA conformation due to elastic and electrostatic forces 
is generally very complicated, and should be treated within the framework of electro-elasticity theory. But 
the available experimental values of elastic stretch constants and the dependence of persistence lengths 
on the ionic strength of aqueous solution, indicate that, strictly speaking, a molecule of the DNA does not 
behave like an elastic rod, The effect could be due to the backbones of the two DNA strands, which 
contain phosphoric groups carrying negative charge that may cause the strands to repel each other, 
facilitate the separation of the strands, and make for the double-helix being less stable, at low ionic 
strength, @. The overall picture strongly depends on the specific arrangement of the constituent base 
pairs, because the forces keeping a base-pair together depend on the choice of constituent nucleotides, the 
number of the hydrogen bonds involved being different, e.g. 2 for adenine - thymine and 3 for guanine - 
cytosine. It is also important, that the relative positions of bases change from pair to pair. Obviously 
we are very far from the simple elastic rod model. Therefore, there is a need for a semi-microscopical 
theory that could accommodate the elasticity of the DNA and take into account its micro-structure. 

The current approach to the problem relies on the hypothesis that it is possible to separate the dy- 
namical modes of the bases and the sugar-phosphate backbone by considering the coupling between them 
as perturbation. The backbone modes are suggested to be strongly overdamped, whereas the modes due 
to the motion of the bases, or inter-strand modes, are assumed to be less sensitive to the viscosity of 
ambient liquid. The theoretical calculations of paper 0] support the statement, whereas the experiments 
on Raman scattering, fl3|.fl2l|, indicate that the attenuation of the modes is substantial, see also [l4| for 
the theoretical treatment of the inter-strand modes. 

The approaches commonly used to manufacture qualitative models that could give an adequate picture 
of the dynamics of the DNA, are generally based on the concept of lattice ( see paper 0, (lQ |. in 
which a coarse-grained model is considered), that is one considers the molecule of DNA as a regular 
structure similar to a one-dimensional crystal. It is necessary to take into account that: (1) the DNA 
comprises the two strands; (2) the strands are bound together by certain forces determining inter-strand 
motions; (3)the relative position of the strands verify the helical symmetry. But, the DNA is not totally 
symmetrical structure due to the choice of base-pairs, which is generally random. This circumstance re- 
sults in considerable theoretical complications. Important, the helical structure supposes the existence of 
the preferred local system of coordinates for every base-pair, and these coordinate systems should change 
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from one base-pair to another, the neighbouring one. Constructions accommodating the phenomenon 
utilize the concept of gauge field, following its current use in field theory and condensed matter physics. 
In the case of the DNA dynamics, it allows for at least qualitative description of its elastic modes. 

Using the lattice has also the advantage of allowing for tautomeric forms of the base pairs. There are 
of two classes of the DNA bases, purine (adenine, guanine )and pyrimidine (cytosine, thymine). Under 
the ordinary circumstances adenine and cytosine are in amino form, and only rarely in the imino one, 
whereas guanine and thymine prefer the keto form, and rarely the enol one. Inside the base-pairs the 
transformation of the tautomeric forms correspond to the tunneling transitions of protons in the hydrogen 
bonds keeping the bases together. It is alleged that the transitions may result in mutations, [la ]. |16j, 
[13 ) [HI ■ The intriguing question is whether there is an interplay among the elastic properties of the 
DNA and the tautomeric transitions. 



II. THE ELASTIC ROD MODEL 



In this section we shall consider more fully the rod model and its possible modifications. One may 
advance the hypothesis that the conformation of the molecule of DNA could be qualitatively described 
by visualizing the latter as an elastic thread that has the elasticity constant, 7, of such a size that heat 
fluctuations due to surrounding solvent be small on a scale called the persistence length. In order to 
make a rough estimate to the effect, we may assume that the elastic energy of a molecule be given by the 
equation 




in which L is the length of the molecule, 4> is the deviation angle of the vector tangent to its central line. 
Then the condition for the persistence length, L p , reads 

E e lastic > ksT 

If we take L « 500A, T « 300 and 4> ss 1 radian, we obtain^ « 10~ 19 erg ■ cm, that is the value which 
is in reasonable agreement with the experimental data, @, [7j- Thus, on the spatial scale between one 
coil of the double helix and 500 4- IOOOA, that is 10 150 base pairs a molecule of the DNA may preserve 
its straight form and has the appearance of an elastic rod. This model of the DNA allows to calculate, at 
least qualitatively, interwound structures called plectonemic supercoils formed by DNA molecules (first 
found by electron microscopy, [l||). 

Within the framework of this model, the double helix of DNA is characterized by three spatial scales: 
(1) the microscale of order 3.4 A, that is the distance between adjacent base pairs along the chain; (2) the 
mesoscale of order 10 3 A, or several persistent lengths; (3) the macroscale of the size of a molecule of DNA, 
that is up to several [im or more. On the microscale, the molecule of DNA is formed by the base pairs 
of purines and pyramidincs linked by hydrogen bonds between the bases; the whole constitutes a double 
stranded structure. The DNA helix considered on the mesoscale, is suggested to have the properties of 
an elastic rod with the torsional and the bending rigidities of about 10 -19 erg ■ cm, Q. The stretching 
of the rod is assumed to be small compared with the bending and twisting, and in many cases may be 
neglected. On the macroscale, a molecule of DNA is flexible, its rigidity does not influence its shape, and 
one can consider it as usual polymer. The whole picture constitutes the so-called worm-like-chain model, 

i- 

The picture of the DNA considered on the mesoscale can be cast in a quantitative form using the classic 
theory of elastic rod worked out by G. Kirchhoff, [20]. It should be noted that the use of the Kirchhoff 
theory for the needs of the DNA involves certain approximations. First, there is a problem of taking 
into account the finite diameter of the molecule, for strictly speaking the Kirchhoff equations are written 
down for the elastic line. Second, the use of continuum mechanics for objects with spatial scales of several 
tens of A may raise some doubts; in fact, there is no continuous medium at hand. Third, Kirchhoff's 
model does not allow for the possible extension of a molecule of the DNA, and the internal degrees of 
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freedom, for example the relative motion of the the strands. Nonetheless, the cautious employment of 
the Kirchhoff theory gives reasonable qualitative results and appears to be justifiable in certain regimes, 
1. 

The static of a system of this kind can be described by means of the effective energy, which constitutes 
the core of Kirchhoff 's model. It is given by the equation 

F = I ds - 2J aijWiUJj (1) 

ij 

in which L being the total length of rod, s the length parameter, its elastic moduli, and uji coordinates 
of a vector that describes the strain of the rod corresponding to the molecule. The vector Q is constructed 
as follows. Consider a local frame defined at a point of the central line of the rod; its first vector v\ being 
the unit tangent vector at the point, and the second and the third, V2, v 3 unit vectors along the principle 
directions of the strain. The three vectors are considered as columns of the rotation matrix R = R(s) 
describing the change of the local frame. The matrix 

w = R~ 1 4-R (2) 
as 

can be visualized as an angular velocity of the local frame, the length parameter s playing the part of 
time, or in the vector form 

w = ^2 (P)jk = -tijk (3) 

i 

The minimization equations for the Kirchhoff energy have the same form as the equations of motion for 
the top (the so-called Kirchhoff analogy) 

d _ _ _ _ dF 

The rich analytical and topological structure of solutions to the top, or Kirchhoff equations has provided 
the necessary framework for their application to conformational problems of the DNA, and the most 
important one has been the theory of supercoiling, [4j. The Kirchhoff model allows for the twist-bend 
coupling, JJ, which corresponds to the chiral character of the DNA. In paper [2l[ the model is further 
extended by taking into account the stretch of a molecule of DNA, and considering the bend-stretch 
coupling. In the notations of papers [U [HI, [HI, (24[ , this means that the elastic energy of the molecule 
reads 



1^ 

E = - J ds [A' ilj + AVl\ + C(n 3 - uj ) 2 + Bula 2 



+ 2DLOo(n 3 -u Q )a + 2G(n 3 -uj )n 2 + 2KujQn 2 a\ (4) 

where the vector Sli , SI2 , ^3 defined by Eqs.([2] - [3]) determines the bend and twist of the molecule, the 
constant Uq, in fact the vector w\ = 0, u>2 = 0, uj 3 — Wq, accommodates the helix twist, and a accounts 
for the stretch. Thus, C is the twist-bending coupling constant of paper [f|, and K the bend-stretch 
one. The conformations of the rod is finally determined by the minimum of the functional given by the 
equation 

M = t~7= - fZ — 2irrLk (5) 
Kb -I 

in which / is an applied tension and Lk is the topological invariant, the so-called linking number, that 
defines the linking between a closed path and its image obtained by a small translation in space without 
self-intersections. The Lagrange multiplier r serves a kind of chemical potential for the linking. In paper 
[22| the authors put forward arguments that there are twist-stretch terms in the elastic energy if the 
molecule of DNA is modelled on a stack of thin, rigid plates that are not permitted to deform. They 
obtain an equation for the elastic energy of the form given by Eq.Q. The problem is related to that of 
small fluctuations of bend in the helical backbone that are important for understanding the mechanism of 
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the torsional stress accompanying the transcription. As is shown in paper (23|, the drag could correspond 
to a torque of 19 -13 dyn cm, which rather surpasses the actual torque involved during transcription. In 
contrast, if the bend fluctuations are not taken into account, the torque is negligible, which contradicts 
the experimental facts. These results are in qualitative agreement with the value 4.5 x 10 -19 for the twist 
rigidity (see also [11]). 

In spite of the successes of the elastic rod model, there has been a considerable criticism levelled at it. 
Baumann et al, @ , have studied the elastic properties of the DNA as a function of ionic strength and in 
the presence of multivalent cations. They measured the extension of the DNA caused by an external force 
applied to it. It should be noted that there are three regimes in the elastic response of DNA molecules, 

1. 0.01 -j- lOpN the molecule behaves as an entropic spring, the worm-like chain model (WLC), [27l |: 

2. 10 -r 65pN deviations from the WLC, enthalpic elasticity; 

3. at about 65pN the molecule suddenly yields in a highly cooperative fashion and overstretches w 1.7 
times, [3(|. 

In the region where the WLC model is valid, a molecule of the DNA is a kind of hybrid of a rigid 
rod and a flexible coil, and is usually visualized as a homogeneous elastic rod. But the assumption is 
in contradiction with the elasticity theory according to which the persistence length P and the stretch 
modulus, S, given by the equation S = EA, in which E is the Young modulus and A is the cross-sectional 
area of the rod, should vary in the same way with the ionic strength whereas the results of Q indicate 
that they change in opposite directions. Another discrepancy comes while considering the Poisson ratio 
a defined by the equation 

B/C = l + <7 

in which B and C are the bending and the torsional rigidity, respectively. The thermodynamical stability 
requires, [2(|, — 1 < a < 1/2, while a < corresponds to the thickening of the rod as it is stretched. The 
values of B are alleged to be 2 ■ 10 -19 erg ■ cm, and C in the range 2 ■ 10 -19 to 3.4 • 10 -19 , [6]. These values 
correspond to —0.4 < a < 0., and therefore the DNA rod should thicken while it is stretched. Thus, we 
should either accept that the DNA is an elastic rod of quite unusual nature, or admit that there is a 
need for a model that could accommodate its double-stranded structure, helical symmetry, and internal 
stacked base pairs. 

Besides the drawbacks mentioned above the rod model fails to accommodate the process of denaturation 
when a molecule of the DNA splits up into separate strands. But, the breaking of the hydrogen bonds 
between the base pairs and the formation of bubbles comprising segments in which the two strands are 
separated, may happen under other circumstances. Equally important, there may exists relative motions 
of the strands which do not result in breaking the hydrogen bonds. In fact, the dynamics of these inter- 
strand modes is accessible to experimental studying, (see for example [l|)j 0, [3l| . [Hj]. It may tell a 
lot about the physics of the DNA. 

One may try to mimic the partition of a single double-stranded molecule into two strands within the 
framework of the rod model, by employing an additional quantity q that indicates a relative displacement 
of the strands from the equilibrium conformation. The procedure is similar to that discussed above in 
connection with the coupling between the stretch and the twist-bend modes of the DNA. What's more, we 
shall see that thermal fluctuations of the field q, result in a contribution to the twist-bend momentum 
and thus turn out to be similar to the stretch term introduced in [2l[, [22[. 

The field q(s) gives displacements of points coinciding in the initial equilibrium conformation of the 
rod, and possible going apart because of the deformation and breaking on the microscopical scale, of 
bonds between base pairs of the molecule. On the mesoscale, which we use to describe the molecule, the 
vectors q(s) are determined at points of the rod by the parameter of arc length, s. Thus, outside the 
region of states where the breaking up takes place, we visualize the molecule of DNA as comprised of two 
clastic wires attached to each other; the whole being an elastic rod that one can bend, twist, and stretch. 
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It is worth noting that on allowing the partition of the strands we at the same admit that the molecule 
can be stretched. In fact, the sugar-phosphate backbones of the strands are formed by strong covalent 
bonds and are hard to be stretched, so that the partitioning of the strands should result in changing 
the total length of the total molecular. Whether it will be diminishing or increasing, depends on the 
deformation of helix that should accompany it. 

The static of a system of this kind can be described within the framework of Kirchhoff 's model of the 
elastic line with the help of its extension with the vector q. It is determined by the energy given by the 
equation, [33[, [3^ |. 

F = { ds I \ ^ ajjUjUJj + huij + \ a [d s g + lu x q\ 2 + U(q) ) (6) 

in which L is the total length of rod, are its elastic modulii, and uji coordinates of the vector of 
deformations described above for the Kirchhoff model. The potential U(q) accommodates forces that 
keep the two strands hanging together. The vector (61,62,63) is analogous to iv of paper [2l| . and 
describes the molecule's winding determined by certain external conditions, for example, a histone. . 
We have used the vector q in the fourth term of equation (J6j> so as to accommodate deformation and 
breaking of bonds between base pairs, and as well as the separation of strands. But here, again, it is 
worthwhile to note that the term is similar to the fourth term in Eq.((4]), see paper [2l[, Eq.(3); besides 
the mathematical similarity we should like to draw attention to the fact that the reason for this lies in 
the relation of q to stretching the molecule. The third term in equation ^ is the covariant derivative 

V s <f = d s q + uj x q 

The covariant derivative is in order because we have to consider the displacements of the strands with 
respect to the local coordinates determined by the helical structure and the strain described by the vector 
q. The procedure is quite common in the theory of gauge fields. 

The state of equilibrium corresponds to the minimum of F and gives the equations 

dU _ 
V^g= V s (m + AqxV s q) = 

oq 

in which V s is the covariant derivative 

V s X = d s X + u) x X 

and m reads 

3 

rrii — 2J a-ijujj 

.7 = 1 



Now let us neglect nonlinear, that is greater than second order terms, in the potential U, and consider 
the effect of thermal fluctuations of the field <f(s), that is inter-strand motion, or in the context of paper 
pH | the stretching, on the total configuration of the molecule. To that end we need to average the field 
q out, and find the effective energy 

e -P F eff = J DqerP F , 0=l/kT (7) 

We aim at a specific, but very important, configuration in which the vector to is equal to the constant 

uj Q = -aT 1 !) (8) 

that provides the absolute minimum for the density of Kirchhoff 's energy given by Eq.(TT]), and the 
conformation of the regular coil for the molecule. Here a is the matrix of Kirchhoff 's modulii a^ . 
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To evaluate the functional integral in Eq.0, we shall employ Feynman's variational principal, j3a |. 
and to that end resolve the expression for the energy ^ in the form F = F + Fx with F given by the 
equation 



F = J ds ( \ £ l? aijUiUij + h 



1 2 



ds 



+ hT, l [B + A[u 2 -^)q 2 



and F\ by the equation 



F x = A J ds I tjjkUiqi ■•^j-qk ^2~,~,'l.'lj 

\ *J* i¥=3 



According to Feyman's variational principal, [35], there is the estimate for free energy T < Tq + (-fi)o 
with T — Ffi uc t of Eq.([7]), and the other averages being 



Dqe 



-0F O 



and 



(Fi)o 



J DqF ie -P F ° 
jDqe-P F o 



The average (Fi)o gives zero contribution owing to the Gaussian nature of the integration. Hence, within 
the limits of accuracy provided by Feynman's principle, we have Ffi uct = T and 



Fef/ = 



+ ^2 bi^i - kT In 



sh(zi) 



with Zi being given by 



1 B 

= - L \l-r + u>$ - wg,., i=l,2,3 



Here B/A evaluates the coupling between the strands of a molecule, and u>oi are coordinates of the vector 
given by Eq.([5]) determining the configuration without fluctuation corrections. 

Let us consider the small coiling and coupling of strands, that is u> <C 1 and B/A <C 1, or Zi -C 1, i = 
1, 2, 3. Assuming uj to be constant, we may cast the equation for the effective energy in the form 



a ij ~ ^TLSij 



UJiiOj 



5> 



24 A 



By minimizing F e ff with respect tow, we obtain the correction of the value of wq effecting the equilibrium 
conformation. Since the contribution of fluctuations is assumed to be small, we may write down the 
corrections to uJq given by Eq. © 



-.fluct 



and 



Sujc] 



--kTLa- 2 -b 



(9) 



From the last equation, we infer that fluctuations of the field q make for an increase in the coiling 
of a molecule. It is easy to estimate a spatial scale on which the arguments given above are valid; 
the fluctuation energy should be smaller than the elastic one, that is counting by orders of magnitudes 
Laoj 2 3> L 2 kTu> 2 which amounts to a/L 3> kT. For the rigidity of order 10 -19 erg cm and room 
temperatures, one obtains the scale of order 5 10 -6 cm, that is the persistence length. To understand 
the estimates given above, let us notice that in the opposite regime, far from equilibrium, we may set 
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B 3> 1, and B/A ^> w 2 . It is easy to see that in this case there are no corrections of the values of C3, and 
consequently no additional increase in the supercoiling of molecule. 



Working within the framework of paper 21|, Moroz and Nelson found the renormalization of twist 
stiffness by bend fluctuations; they averaged the partition function corresponding to the functional given 
by Eq.Q, obtained the Schrodinger-like equation for the correlator of the tangent vector to the molecular 
axis, [24| . [25j . and found the torque, r(/, a), as a function of applied stress / and external twist er = ALk 

T{f,0) = (10) 

C- 1 + UAyjAf/k B T 



Since the setting of the two problems is not identical, we are not in a position to compare Eqs.© and 
p0|) : qualitatively, they are related to the same phenomenon of renormalization of stiffness by bend 
fluctuations. 

It is worth noticing that the hydrogen bonds of DNA can be deformed and even broken due to the 
local action of an external agent, for example enzyme. We may try to describe the process by an energy 
term that be included in Eq.(fT|) 

L 

Fir* = / ds (0 ■ q + yyq) 5 (s - s ) (11) 



with [] and 7 being a constant vector and a matrix, respectively. The <5-function factor is used for 
describing point-like action of the external disturbance. Consider the case of the initial value for b equal 
to zero before the breaking up of a molecule. For small q, using the linear approximation, we have the 
following equation for q 

d 2 (3 1 

— s q-n 2 q= — S (s - s ) + —•yqS (s - s a ) (12) 
with fl 2 = B/A . Up to the second order terms the minimization equation for lu reads 

a VcD - A V (jj^q* <t\ = 

From the last equation one can infer that the external action due to Eq. (|ll[) results in the formation of 
an effective rotating moment b ext that reads 

bext = ~^— L -r-G (s, s ) [J3 x 7/3J 

with G (s, sq) being the Green function of Eq.f [T^|) . 

From the equation for G(s, Sq) it is easy to see that the size of the defect changes from L, or the 
mesoscopic scale, for small O, to I/O for large fi. 



G(s,s ) = - _ . , n = . (9(so - s)sh(Qs)sh ((L ~ So )Sl) 

+ 6(s- s )sh(Qs )sh({L- s)Q)) 

with P and 7 being a constant vector and a matrix, respectively. 

The local influence due to a chemical-biological agent may result in the strands of a molecule being split 
up on a scale of order 1/Q determined by the potential describing the hydrogen bonds, so there is a kind 
of Lindemann's criterion: the breaking up of a molecule takes place if y/ B /A cx 1/L. The magnitude of 
the splitting is determined by the properties of the agent. A drawback of the model which is based on a 
quadratic Lagrangian, is that it does not allow for a threshold effect, which might be taken into account 
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by considering high order terms. In the region of conformations far from the splitting of a molecule 
the fluctuations of the field q(s) are suppressed by the bond due to U(q), and there is no enhancing 
the supercoiling, in contrast to the region close to the splitting, where the potential U(q) is effectively 
small. Since the minimization of energy makes sense only on the mesoscale of several persistence lengths, 
the conformation of a molecule being determined by entropy on the macroscalc , different segments of 
a molecule should be broken up independently from each other so as to give a blurred character to the 
transition corresponding to the splitting up of the molecule. 

There are various arguments to the effect that the double-stranded DNA should suffer local denatu- 
ration and open up locally so that hydrogen bonds between base pairs be broken. Besides physiological 
processes that involve the 'unzipping', it can take place spontaneously due to fluctuations owing to the 
small energy required, less than 3fcsT, Q. Therefore, the breathing of DNA is a phenomena that could 
happen in many situations. 

It was G.S. Manning, [36|, who suggested that the real flexibility in the DNA molecule could be due 
to opening of base pairs and the breathing fluctuations be related to elastic properties of the DNA. The 
solution to the problem requires the knowledge of the rate of the breathing fluctuation, and at this point 
there is some controversy. 

By now there are two methods for studying the breathing fluctuations: the NMR and the fluorescence 
correlation spectroscopy. 

The NMR measures the exch ang e of protons from imino groups with water, which are suggested only 
to occur from open base pairs, |37l ] According to paper [38j the opening of base pairs may require prior 
unwinding or bending of the DNA double helix, which does not necessarily lead to imino proton exchange. 
The life times of base pairs and open states depend on temperature and bases involved. In the B-DNA 
at 15°C typical lifetimes range 0.5 through 7 ms for A ■ T pairs and 7 through 40 ms for G ■ C . The open 
state life time is 10— 100ns, [38|. Leroy et al, [38|, estimate the activation enthalpy for C4-G5 opening to 
within 45-^56^7/7710?, or 10.7-j- 13.3Kcal/mol . It is important that the NMR of imino proton exchange 
measures the lifetime of a single base pair, (39j . 

The fluorescence spectroscopy relies on the translation of base pair fluctuations in fluorescence fluc- 
tuations. Altan-Bonnet et al, [4(|, use synthetic DNA samples containing modified bases tagged with a 
fluorophore and a quencher. When the DNA structure is closed, the fluorophore and the quencher are 
in close proximity and the fluorescence is quenched, it is again restored if the structure opens so that 
the fluorophore and the quencher are pulled apart. Thus, the fluorescence spectroscopy detects the local 
denaturation of the DNA, or bubbles of 2 to 10 base pairs with lifetimes in the 50/is range at 37°C. 

The divergence between the NMR and the fluorescence spectroscopy measurements could be explained 
by several reasons. First, the NMR picks up very fast modes owing to the imino proton exchange being 
very sensitive to the conformation of a base pair, whereas the fluorescence spectroscopy picks up large 
scale deformations of the DNA. Second, the NMR measures the life time of a formed base pair; the 
fluoresceence spectroscopy measures the life time of the open state. Third, the range of temperatures and 
the composition of the strands in the two experimental settings is different. For these arguments I am 
indebted to G. Altan-Bonnet. 

Thus, the available experimental data indicates that partial denaturation of the DNA is a common 
phenomenon, and it should be taken into account while using the elastic rod model. The latter requires, 
generally, serious modifications when the internal motions of the double helix need accommodating. 

There are arguments of qualitative nature in favour of the fact that the fluctuations of open hydrogen 
bonds result in an effective interaction between the strands of the double helix (cf. p(3] about the bonds 
between base pairs and those with water). We shall describe a hydrogen bond as two- level system, |0) 
and |1), with the energy difference e = Eq — E\, and the Hamiltonian eo"3 where 03 is the third Pauli 
matrix 



1 

0, = - 



1 

-1 



To describe the coupling between the Pauli operators of hydrogen bonds between base pairs and the 
classical system due to the elasticity of the molecule, we shall employ the field q (see pj6j). We suggest that 
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the interaction between the elastic forces and a single hydrogen bond is small, and therefore perturbation 
theory could be employed, that is a kind of Born-Oppenheimer approximation in the sense that one part 
of the system is considered to be classical, whereas the other one quantum. We cast the equation for 
energy in the form 

E = F(iU,q)+E mt (q,a)+E exc (13) 
where F is the energy of elastic rod given by Eq.([5]), and E exc is the energy of hydrogen bonds, 

AT-l N-i 



n=0 



with <r„ being the matrices 



0"n = \{^n ±i(J l) 



The second sum in the above equation accommodates the possible propagation of opened base pairs. The 
third term in Eq. (|13[) is the interaction between the classical elastic system and the quantum one given 
by the hydrogen bonds. The field q has the sense of a mean field describing the partition of the strands 
and its characteristic spatial scale is hundreds of A. The spatial scale of the Pauli operators is a few A. 

It is important that the operators er^, er+, cr~ are intimately related to the conformation of the molecule. 

In fact, the breakdown of a hydrogen bond follows a certain direction in space. At each site n cor- 
responding to a base pair of the molecule, there is a local frame formed by unit orthogonal vectors 
wi, u>2 W3, for which the vector wi is tangent to the axis of the double helix, W2 is normal to the axis, and 
W3 indicates the direction of the bond's breakdown. The operators cr^, <r+ , a~ indicated above are chosen 
in accord with the frame w±, W2W3. Thus, the term cr^a~ +1 in the energy E exc involves the operators 
that act in spaces (|0) n , |l) n ) an d (|0)„+i, |l)„ + i). To form the interaction term we need to cast q and 
the Pauli operators in the same co-ordinate system. By using the unitary transformation 

we may find the operators s l n , and specifically s^, that correspond to a l n in the laboratory co-ordinate 
system and describe the state of the hydrogen bond at site n from the point of view of an external 
observer. In fact, there is the equation 



U^alU = Rika* 



fc=l 



in which the matrix Rik corresponds to the rotation that brings the frame wi, W2W3 in the standard 
laboratory one. Therefore, we choose the interaction energy in the "minimal" form 



7V-1 



E in t(q, 3) = 7 X! 5" ' Rn °" ( 14 ) 

n=0 

We shall confine ourself to the circular conformation of the molecule for which the matrices R n have the 
form 




Therefore the interaction terms read 

q n ■ Ra n = (cos (j> n ql + sin<£ n g£) a \ + (- sm<p n ql + cos<j> n q%) 

and we have 

iV-l 

Ei + E exc = 7 J2 [(cos 4> n Qn + sin (j> n q%) a 2 n + (- sin $ n q\ + cos <f> n ql)] 

n=Q < +PJ2n=0 [ a n <?n+l + °n <+l) 
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Since we are considering the molecule of DNA at temperature far from denaturation, and therefore in 
accord with the results of papers [3?], HE the number of excitations, that is of broken hydrogen bonds, 
is small. Consequently, we may utilize the method of approximate secondary quantization, [281 ] , familiar 
in the theory of spin systems. Let us consider the Bose operators bk, b^ 

[h, b m ] = 0, [b~£, 6+] = 0, [b' k 6+] = 8 km 

By using the substitution 

<rl = 2- b k b k, crfe=6^, ,a k h b k = b k 
we may cast the energy Ei nt + E exc in the form 



eN 

Ei + E exc = — - + -^] (-sin cj) n ql + cos <j> n q^) 

n=0 

JV-1 JV-1 

+ [ _ 7 (- sin <f> n ql + cos cf> n ql) + e] b+b n + (3 ^ + b n b n+x) 

n=0 71=0 
. JV-1 

+ ( cos< ^ + sm<j) n ql) (b+ - b n ) 

n=0 

The first line in the above equation corresponds to the 'vacuum' fluctuations, which may be of interest 
as regards the conformation of the double helix given by the matrices R n and the field q. The above 
equation is quadratic with respect to the Pauli operators bk , b^ and we can cast it into the diagonal form 
by employing Fourier transform 

f_i 2iirnk 



B, - £. N h 

v n— 

N _l —2innk 
B+ = -L V 



so that the energy -Ei n t + -E ea; c be given by the equation 

JV-1 

+ £ CTC = const + - ^ (- sin 0„ + cos 0„ g 3 

n=0 

AT-l JV-1 

+ E ^ lk2 B+B k2 + (A n B++A* n B n ) 

k\k2 71=0 

in which and A„ read 



e kl k 2 = 4 lfc2 (e + 2/3cos— (15) 



2i7r(fc2 — fei)n 

■1 ^ AT - 1 (- sin ^g* + cos <M 3 ) e 

71 = 



4 - I^2_V 





JV 



2 VF £: 



-l -2i7r(fc 2 - fci)n , 2 . 3 s 

p ^ (cos <j> n q n + sm<j)nq„) 



We may eliminate the linear terms in B n , B+ using the canonical transformation 

B n — ► S 1 B n S, B% — > S 1 1 5 
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or more explicitly 

B n — ► C„ = B n + l n , £>+ — > C n = £>„ + i* 

We obtain the equation 

N-l 

Ej + E exc = const + E con f + ^ e k 1 k 2 C£ i Ck 2 

k\k 2 

in which the conformation energy E con f reads 

N N-l 

E conf =2^,-1 (-sin<i> n ql + cos <f> n ql)+ 3 (e _1 ) fclfc2 A lA*k 2 (16) 



According to the results of Leroy et al [37|, |38|, |39| the effective temperature of the excitations due to the 
opening of hydrogen bonds is very low, and therefore the energy of the system is given by E con f. As was 
assumed the interaction between the elastic part of the total energy and that due to the breakdown of the 
hydrogen bonds ia also small, so that we may use the approximation of Born-Oppenheimer. Therefore, 
it is possible to consider the operators B n , B+ assuming that q and uj are constant. At the same we 
shall assume that 7 is small and make all calculations up to the order j 2 . In fact, we are considering the 
interaction energy given by the term 7 q n ■ Ra n , which is smaller than e. Since A n cx 7, we shall neglect 
terms of order 7 2 in {e )kik 2 an d obtain 

I -\\ $kik 2 

^ > k 1 k 2 ~ 2nk 
e - 2/3 cos — 

On setting q n = q 2 + i q^ we cast the conformational energy in the form 

N-i o 2 N ~ 1 1 i^n m2 



where 



2nk 

n=o ' A " n=o a — 2/3cos 

N 



v 1 2mmk 
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D(k) = — F Ve N \cos< 



m=0 

On integrating out the phases <j> n , we find the effective potential 



3e 2 

U(q) = E conf = — ^ \ln\ 



32e 

n=0 



or using the integral form 

TT(n) = E c = 

32ed 



U(q) = E conf =^- I \q n \ 2 ds (17) 



Within the framework of the model considered above the excitation of the hydrogen bonds, i.e their 
breakdown, is visualized as an ideal gas, similar to that of magnons. It is applied only to regimes far 
enough from the denaturation, when the number of the broken hydrogen bonds is small, [33] • Nonethe- 
less, it sheds some light on the formation of locally denaturated regions, or bubbles, of the DNA. It is 
worthwhile to recall that G. Manning had suggested that there are two types of excitations of the DNA; 
the modes of bending and those of breathing, (36|. By now it is generally accepted that the bending 
modes are of high frequency, and could be related to nonlinear phenomena in the DNA. The breathing 
modes according to [37| are in GHz-rcgion, which is at the edge of the elastic modes of the DNA. If we 
direct our attention primarily to bubbles formed by opened base pairs, characteristic time appear to be 
by orders of magnitude larger, ■ It is also worthwhile to note that the model of elastic chain of the 
Pauli operators, has some bearing on the opening of base pairs and the breakdown of hydrogen bonds, 
mentioned on pJ3J 
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III. THE EXTENDED LATTICE MODEL 



As was discussed above, the elastic rod model is not generally sufficient for describing the conforma- 
tional dynamics of the DNA. Important, it does not allow for its intrinsic degrees of freedom corresponding 
to the structure of the double helix. To build an adequate model to the effect, is a difficult problem, 
and an attempt to manufacture it 'from first principles' is doomed to failure. Thus, there is a need for 
drastic simplifications, and it is necessary to take into account: (1) the DNA having the two strands; (2) 
the base-pairs being linked by the hydrogen bonds; (3) the helical symmetry of the DNA. The problem 
still waits its general solution, but specific cases are nonetheless tractable. In this section we are going 
to see what could happen if elastic modes that can be expected within the rod model, may interact with 
internal motions of the double helix. 

We consider short pieces of the DNA, of several persistence lengths, so that the spatial conformation 
of the molecule on the mesoscale , is not of primary importance. We focus on the internal dynamics, 
trying to accommodate the above requirements through a one-dimensional lattice model of the DNA. The 
key point in this respect is the wise choice of dynamical variables that could give a picture of the DNA 
dynamics, both simple and adequate. El Hasan and Calladine, 1], give the framework for such analysis 
by setting up the scheme for the internal geometry of the double helix of the DNA. They describe the 
relative position of one base with respect to the other in a Watson-Crick base-pair, and also the positions 
of two base-pairs, by introducing local frames for the bases and the base-pairs, and translation-slides 
along their long axes. 

We follow the guidelines of paper [l[, but aiming at a qualitative description of the DNA dynamics 
use a simplified set of variables. We shall describe the relative position of the bases of a base-pair by 
means of the vector Y directed along the axis of orientation for complimentary bases inside the base pair; 
Y being equal to zero when the base-pair is at equilibrium. The relative position of the base-pairs is 
described by the torsional angles <j> n , which give deviations from the standard equilibrium twist of the 
double helix. Thus a twist of the DNA molecule, which does not involve inter-strand motion or mutual 
displacements of the bases inside the pairs, is determined by the torsional angles 4> n that are the angles 
of rotation of the base-pairs about the axis of the double-helix. The twist energy of the molecule is given 
by the equation 



in which / is the moment of inertia, and r is the twist coefficient, which for the sake of simplicity and 
taking into account the qualitative picture at which we aim, are assumed the same for all the base-pairs. 
Inter-strand motions should correspond to the relative motion of the bases inside the base-pairs, therefore 
the kinetic energy due to this degree of freedom may be cast in the form 



where M is the effective mass of a couple. 

For each base-pair we have the reference frame in which z-axis corresponds to the axis of the double 
helix, y-axis to the long axis of the base-pair, x-axis perpendicular to z- and y- axes. At equilibrium the 
change in position of adjacent base-pairs is determined only by the twist angle of the double helix. We 
shall assume f2 = 27r/10. To determine the energy due to the inter-strand displacements we need to find 
the strain taking into account the constraint imposed by the helical structure of our system. For this end 
one may utilize the method employed by G.Kirchhoff for the twisted rod, that is the covariant derivative, 
as was done in paper [4l( for the DNA molecule. But a more simple and straightforward approach is 
possible. 

Let us confine ourself only to the torsional degrees of freedom of the double lattice and assume the vec- 
tors Y n being parallel to x-y plane, or two-dimensional. Consider the displacements Y n , Y n+ \ determined 
within the frames of the two consecutive base-pairs, n, n+1. Since we must compare the two vectors in 
the same frame, we shall rotate the vector Y n+ \ to the frame of the n-th base pair, 
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n 



Y^ k = R- 1 {<t>)Y n+1 
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Here R 1 (4>) is the inverse matrix of the rotation of the n-th frame to the (n+l)-one given by the equation 



R(<f>) 



cos <p — sin < 
sin <b cos 6 



(18) 



The matrix R is 2 by 2 since the vectors Y n are effectively two-dimensional. Then the strain caused by 
the displacements of the base-pairs is determined by the difference 

\r back \} 

I n+l x n 

For this argument I am indebted to D.I. Tchertov. 

It is important that the angle <f> is given by the twist angle, £1, describing the double helix, in conjunction 
with the torsional angles <£>„, so that 

4> = fi + 4> n+ i - <j) n 

Therefore, the energy due to the inter-strand stress reads 

K 



E 



M 4 2 



R l (Q + <f> n+1 - 4> n )Y n+1 -Y„ 



It corresponds with the fact that the equilibrium position of the double helix is the twisted one deter- 
mined by O and all <j> n being equal to zero. We suppose that the size of DNA molecule is small enough 
that it can be visualized as a straight double helix, that is not larger than the persistence length. Hence 
the number of base-pairs, N < 150, approximately. Combining the formulas given above we may write 
down the total energy of the DNA molecule in the form, [42J , 



H = 



E 

n 

E 



2 



^2 _ ^«) 2 



K 
2^2 



R-\n + 4> n+1 - 4> n ) Y n+1 - Y n 



2 " 



(19) 



in which K and a are the torsional elastic constant and the inter-pairs distance, correspondingly. In 
summations given above n is the number of a site corresponding to the n-th base-pair, and n = 1, 2, . . . , N, 
N being the number of pairs in the segment of the DNA under consideration. The last term, e/2Y 2 
accommodates the energy of the inter-strand separation due to the slides of the bases inside the base- 
pairs. 

It should be noted that the dynamical variables <j) n and Y n are of the same order of magnitude, that is 
the first. Consequently, preserving only terms up to the third order, we may transform Eq. (|19p . so that 
it takes on the form 



H 



E 



2a 2 



fn+l 



2^ 1 9 Yn + 2a 2 



J R- 1 (fi)y„ +1 -y„ 



-r 2 

2 n 



— E 



) R-\n)Y n+1 x Y n 



We have used the fact that the axis of the double-helix is directed along Oz-axis. 

Let us simplify Eq. (|20[) by diagonalizing it with the help of the unitary transformation 



Y n S u n , 



S 



V2 V2 
i 1 

V2 V2 



(20) 
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which is a two by two matrix, for the vectors Y n and u n are effectively two-dimensional, their third 
coordinates being equal to zero. The equation for the energy ([20]) takes on the form 



* = E 

n 

+ E 



1 12 T 

2^ + 2^ 
M u u * e _ _ AT 



in 2 _ 2 

"n+l "n 



A" 



E 1 



fn+1 - <Pn 



1 2 

-ie in u] l+1 u n +ie~ m u 2 n+1 u n 



The star * signifies complex conjugation. 

We can further simplify the equation for the energy by applying the Fourier transform given by the 
equations 



n=+N/2 

f*=-m E ei " 0fl ^ 9= ^ m; m = > ±l J ±2,...,±-; 



n=-N/2 

It is important that after the Fourier transform the variables u n verify the following equations for their 
complex conjugates 

% = iu ~qi "<j= ( 21 ) 

The equation for the energy can be cast in the form 



n 



E 

q 

E 

q 

K 



1 X X* , r • 2 «9 

— <Pa <t>„ H n sin — 

2 Vq Vq 2a 2 2 



'q v q 



M ■_ ■_* e _ _ 2K ( 2 n-aq , »i 
yU g • u g + -u g • u 9 + — | sm •■ " 



g'g" 



u g M g + sin 



2 tt + aq 2 *2 



q> "-q'+q" 



(22) 



in which 

<7* = Q/a 



The above equation serves as well a Hamiltonian that describes the dynamics of a molecule of the 
DNA, within the framework of the present model, up to terms of the third order. It is worth noting 
that the latter is based on the assumptions given above concerning the basic structural properties of the 
DNA, that is it takes into account its two-strand structure, the lattice formed by the base pairs, and the 
helical symmetry. The specific feature of the Hamiltonian is the presence of the interaction term that 
describes three- wave interaction, [42] , and may result in resonance. We shall utilize the fact for deriving 
the parametric maintenance of the u q modes, i.e. the HBS (hydrogen-bond-stretch) modes, (see below). 

One can obtain, in the usual way, the equations of motion for , a = 1, 2 and <p q , from the equation for 
the energy indicated above. The essential point is the effects of dissipation, which are due to ions in the 
close neighborhood of the molecule and water effects, see [43| . The dissipation could be accommodated 
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by writing down terms linear in ii®, <f> q . We shall take into account external force, or torque T q , only in 
the equation for (f> q , for it corresponds to external degrees of freedom of our model. Thus, the equations 
of motion can be cast in the form 



u 



« + ul qUq * + lu u" + ^ W ^><- q > = 0, (23) 



Here 



2 4K . 2 n+{-l) a aq e 2 4r . , aq 
<=M*« m 2 + M> ^ = 7^ Sm T (25) 

are the dispersion laws for the fields u q , a = 1,2, and <p q . We see that the spectrum of 4> q has a typical 
acoustic character, whereas that for u q has a local minimum determined by the helical twist, Q. Thus, 
the spectrum of our model is in qualitative agreement with conclusions of [14j . The specific nature of the 
torque is to be specified elsewhere, for the moment, we shall consider general dynamical phenomena to 
which the torque may be conducive. 

Let us suppose that for one thing the amplitudes of the HBS-modes given by u q be so small that the 
quadratic term in Eq. (|24[) can be neglected, and for another the external torque T q be appreciable enough 
to maintain the vibration of the torsional mode <p q . Thus, we can visualize the latter as a pump mode 
that interacts with the HBS-mode u q through the non-linearity in Eq. (|23]) . We shall confine ourself to 
the case of the torque T q being non-zero only at q = and having a frequency 2uj. Therefore, the forced 
wave, or the pump wave for the HBS-mode, has the form 

= $ Sgg, , </>_f = e~ a ^ $* S-gg* (26) 

To obtain larger values for the pump wave, <f> q , the resonance condition 

LU q * = 2lo 

should be verified, even though the resonance behavior of the torsional (f> q -mode itself could be attenuated 
by dissipation, i.e. it may be a mode of small amplitude. 

The equations of motion for u q in the pumping regime read 



i) a I .-.: 2 Ti a I -.- ii a I ^ K (Ar l2ujt n a + A* p- i2ujt n a \ — n 

U q +UJ aq U q+luU q +J^- 2 -J={Ae U g _ ? „ + A e U q+q *) - U 

here 

A = er in $ 

Note that the momentum conservation in the q-values is preserved, as required by the three-wave inter- 
action. The equations indicated above can be cast in the matrix form 

'U a + Q? a u a + j u ii a = (e !2wt K + e-' 2wt K+) u a (27) 

where K, and tC + arc hcrmitian conjugate, and 

/ 2K smQ\ 2 , , l9 
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It is worth noting that Eq. (|27|) is a kind of matrix Mathieu equation. In fact, we can apply to it Rayleigh's 
method for studying parametric resonance, (44[- For this end let us look for the solution to Eq. ({2"T)) in 
the form of a series 



u{t) = A x e Mt + Bi e- luJt + A 3 e lZu3t + B 3 e~ l3ujt 



On substituting the expression given above into Eq. (|27j) and preserving only the terms corresponding to 
e ±luJt : we obtain the equations 

[{-lo 2 + i lu Lo)l + Lo 2 a ] Ax + KB X = Q 
[{-lo 2 -i 1u lo)1 + Lo 2 a ] Bi + /C+Ai=0 

The compatibility condition of the equations indicated above can be cast in the form of determinant 
for the block matrix 



det 



JC+ 



= 



(28) 



Here Co 2 is the matrix of frequencies given by Eq. (|25|) . and lo 2 and "f u u) are the scalar ones. We can 
transform Eq. ([28|) into a more amenable form. Let us notice that it is equivalent to the equation 



det 



U> 2 —UJ 2 + ij u Ld JC 

K+ "< 2 - >■» 



W a — W — IJuiV 



1 -{uj 2 a -Lj 2 + iy u v) 1 
p 2 JC+ 



in which 



Ma z VN 



P = 



2K sinfi 

and the matrices JC + and K, verify the equation 

-1 + p 2 JC JC+ = 

We have used the fact that for the range of frequencies we are considering, the matrix 

LJ 2 a - UJ 2 + i"/ u LU 

is not degenerate. Therefore, the equation given above is equivalent to the following one 

det [-1 + p 2 {ujI -lo 2 + ij u to)J{uj 2 a - uj 2 - i-/ u lu)J+] = 
in which the matrix J is given by 

Jqq< — S q ' q -q t 

We may cast the last equation into the form 



(29) 



,o 2 - \ / 2 2 - \ ( ^ n ^ \ iii'"' 

[U aq -LU -t~f u Lu){LJ aq _ qt -UJ +IJ U UJ) - I J^2—/ = ) \ A \ =° 



(30) 



which is quite similar to the usual condition for parametric resonance. Solutions to Eq. (|30[) are generally 
complex and therefore correspond to attenuated regimes. But there is a specific wave number, q res , for 
which the solution gives the real frequency lo, and it is easy to see that it should satisfy the constraint 



aq-q. 



} aq' 



at 



q = q r 



(31) 



19 



Thus, we may cast the condition for parametric resonance in the familiar form, [44j |. 

<" 2 -<..) 2 ^ 2 " s -(ww) 2|A|2=0 (32) 



The existence of the parametric resonance discussed above presupposes that the bands of the torsional 
acoustic (TA), the (f> q modes of the present paper, and the hydrogen-bond-stretch (HBS) modes, that 
is the Uq modes, of the DNA interpenetrate each other. So far, the knowledge of the spectra of the 



DNA vibrations relies mainly on the computer analysis within the framework of quantum chemistry, 14 1 . 
According to these results, it is accepted that the conformational dynamics of the DNA is confined to 
elastic vibrations of the DNA molecule in the range of 10 9 ±10 12 iiz, g|. Kim and Prohofsky, [3, 

claim 

that the region comprises two domains, which correspond with different degrees of freedom of the molecule: 
(1) acoustic modes, which do not involve the hydrogen bonds; (2) modes that stretch the hydrogen bonds 
between the base-pairs (the HBS modes). Local minimum of the frequency is characteristic of the HBS- 
modes, |14 |; its position depending on the choice of the band. The vibrations of the DNA, which are 
ascribed to the inter-strand modes, were observed in the low- frequency Raman scattering, [l2 |. [l3j . and 
the Fourier-transform infra-red absorption experiments, [45|. Globus et al, [46j], report the existence 
of internal modes generated by the interaction of artificial DNA-type molecules with electromagnetic 
radiation in sub-millimetre range. It should be noted that the type of modes observed depends on the 
kind of DNA samples, i.e. in aqueous solutions, or films and filaments, |46| . The experimental data, 
[ill, is not conclusive as to the relative positions of the acoustic and the HBS modes. The inter-strand 
and the acoustic modes of the DNA are alleged to be overdamped, [l4j. But the opinion of scientific 
community in this respect is not unanimous. First of all, they are observed and measured, fairly well, 
in experiment, |4?|, [l2j], [4!|. Second, according to paper Q it is the modes related to motion 

of the sugar-phosphate backbone that are overdamped, that is the <j> q modes discussed above, whereas 
the inter-strand modes, or the above u q are not. Third, Davis and VanZandt, [43| . had shown that the 
dissipative effects suffered by the modes are greatly diminished in case the GHz-frequency range is taken 
into account; in the re gion it is necessary to employ the so-called Maxwell model of hydrodynamics, or 
the Leontovich theory, [47| . Thus, we see that there are serious arguments against the inter-strand modes 
being overdamped and absent. 

We are in a position to assess the action of mw-radiation on the molecule of the DNA. The question 
which has aroused controversy during the last fifteen years. We wish to make it clear that there is no 
" scare- mongering" , but merely a suggestion to employ microwave radiation for studying the biophysics 
of the DNA. In fact, we feel that the technic of Raman scattering, which had been successfully used for 
detecting the inter-strand modes, in conjunction with the microwave radiation could be instrumental in 
studying intrinsic motions of the DNA. 

The key point of the theoretical analysis of the interaction of electromagnetic radiation with the DNA 
is accommodating the fact that the wavelength of radiation is by many orders of magnitudes larger than 
the characteristic size of the region of the molecule involved in the process. It was Chun-Ting Zhang, 
[5C| . who suggested a mechanism to overcome this difficulty. The main point of Zhang's argument is that 
the helical configuration of the electric dipoles corresponding with the base-pairs makes the interaction 
of the dipole P and the field E 

U = -P-E 



dependent on angle. Therefore, different torsional momenta are applied at the base-pairs. The equation 
for the energy of interaction between the dipoles of DNA and an incident micro-wave reads 

- E ■ R(nn + <j> n )P 

n 

Here R(n£l + cf) n ) is the rotation matrix given by Eq. ([T8|) . and P is the dipole at site n = 0. Consequently, 
even though on the molecular scale the radiation has a plane wave configuration, it still twists the DNA 
molecule about the axis of the double-helix. Since the momenta changes periodically in time with the 
incident wave, the irradiation results in a periodic stress that may produce elastic vibrations in the DNA 
molecule. Zhang suggested that the force may generate resonance vibrations, resulting in a cross-over 
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mechanism which takes up initial torsion excitations and transforms them into longitudinal acoustic 
vibrations. 

In the present paper we will try to combine Zhang's mechanism, (50j , and the excitations of the double- 
helix studied by Prohofsky and Kim, [T3| , with the view of generating inter-strand waves in the DNA by 
mw-irradiation. In contrast to the original idea by Zhang, we do not utilize a cross-over into longitudinal 
acoustic vibrations, but employ the interaction between torsional oscillations and the inter-strand ones, 
i.e. the three- wave, given by Eq. ((22| . 

The main point is that by expanding the rotation matrix i?(nf2 + <j) n ) in the angles <p n and keeping 
only the first order terms, we may cast Zhang's interaction in the form 

Hz = - <t> n (E x P n ) 3 + const, P n = R(nQ) P Q (33) 

n 

in which P is the dipole vector at site n = 0. Next, by using Eq. (fT5]) for the matrix R(nfl) and neglecting 
the constant term we may cast Eq. (|33[) in the form 



Hz = \ E { einU [(^ x - ifi ■ p o)] + e ~ m ° [(^ x p oh + »(^ ■ Po)]} 

n 

On applying the Fourier transform for the <p n , and utilizing the equation 

_L V e i{n±aq)n - 8a 4. 

n 

we obtain the following expression for Zhang's interaction 

N ( - - - - --I 

Wz = j {0 9 =a [(E x P Q ) 3 - i(E ■ P )] + <f> q ^a [(E x P ) 3 + i(E ■ P a )]j 

Hence, the torque T q in Eq. (l24|) corresponding to TLz is given by the equation 



2 r Z 

in which 



T = - 5 q ,- 9 . + — 8 q ,f q, - n/o (34) 



2 



(E x P ) 3 + i{E ■ P ) 



It should be noted that are the local minima of the HBS-modes. From Eq. (j25]) we infer that q verifies 
the constraint given by Eq. (|31[) reads 

Qres = ^ q* (35) 

It is worth noting that the wave numbers and q res correspond to the wavelengths of one and | turns 
of the double-helix. 

Equations given above provide an opportunity for making numerical, order of magnitude, estimates, 
which enable us to assess the effect of mw-radiation on the HBS-modes. From Eq. ([34|) we infer that the 
torque T has the size 

Tote 2iuJt EP 

where E and P are the external field and the dipole moment of the base-pair, respectfully. Next, suppose 
that the resonance condition 
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be true, so that the action of the radiation on the torsional modes should be the largest possible. Then 
the amplitude of the pumping wave, <j> qt , according to Eq. (|24[) , is of the order 



VN EP 

$ oc -V- 7; (36) 

Next, we turn to Rayleigh's condition for the parametric resonance of the HBS-mode given by Eq. 
For the pumping wave corresponding to Eq. ([36|) . it gives 

2 
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Hence we have the threshold 



2^8^^ EP , . 

7 " 7 ^M^— (37) 

which is the condition that the energy supplied to a DNA molecule is greater than that dissipated, so that 
the maintaining of the HBS-mode can take place. We suppose that the frequency of the HBS-modes, as 
given by Eq. (|25|) . is generally determined by the gap term e/M in the equation for u aq and the first factor 
in Eq. (|37[) does not differ much from unity. It signifies that the energies of the inter-strand separation 
per base-pair and the twist of the relative positions of the two adjacent base-pairs, should be comparable. 
At any rate, the hypothesis appears not to contradict the data reproduced in paper |14| . If so, we could 
have the estimate for the dissipative constants, at least by orders of magnitude, 

EP 

lul4,<— (38) 

On utilizing the relation 

E cx 2 

which follows from the expression for Pointing's vector 

S = —E x H, 

in which c is the velocity of light, we cast the estimate given by Eq. (|55|) in the form 

P I S 

7«70<2y^ (39) 

in which S is the power density of the interaction. If we assume 

P cx 1 Debye or 1(T 18 CGS 

and the inertia coefficient / cx 10 -36 grcm 2 , corresponding to the mass of the base-pair cx 10 -22 gr, and 
the size cx 10 A, then for the power density S cx 100 mWf cm 2 , we have 

lu 70 < 10 16 Hz 2 or 7« , 70 < 10 8 Hz 

The estimate suggests that the effect produced by mw-radiation is to be looked for at the edge of the 
GHz zone, for in this case the requirement on the line-width is less stringent. It should be noted that the 
crucial point in assessing the feasibility of experiments on mw-irradiation of the DNA, and its possible 
influence, is the part played by ambient solvent and ions contained in it. In fact, the irradiation may result 
in just heating the solvent, so that the dissipation due to the ions takes up all effects on the molecules 
of DNA. Generally, the thin boundary layer of water and ions close to the DNA-molecule may have an 
important bearing on the dynamics initiated by the incident mw-radiation and result in the overdamping 
of the molecule's torsional oscillations. 
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Davis and VanZandt, 43], put forward arguments that the ions contained in a layer close to the DNA 
molecule should have an influence small enough to allow the survival of the effect due to mw-irradiation. 
The part played by the dissipation caused by water is more subtle. 

The current arguments, [STj] , about the overdamping of the DNA elastic modes, rely on the Stokes law 
for frictional force, F = 6Trr)Rv, for a sphere of radius R moving in a fluid of viscosity rj at speed v. It is 
important that for the specific case of the DNA it should involve the GHz region of frequencies. But, the 
classical hydrodynamics, that is the Navier-Stokes theory, breaks down in the region, as can be inferred 
from the phenomenon of light-scattering in liquids, which is characterized by the triplet structure: the 
central Rayleigh line, v, due to the elastic scattering, and the Mandelstam-Brillouin doublet, v ± /, of 
the inelastic one; / being the frequency of elastic waves in liquid. The classical hydrodynamics gives the 
width of a line in the Mandelstam-Brillouin doublet larger than the distance between this line and the 
maximum of the central line of the triplet, so that the discrete triplet structure should not be observable; 
in fact, it is, [13]. Mandelstam and Leontovich, brought about the solution to this problem by 
using the relaxational theory of hydrodynamics in which liquid is considered as a viscous elastic medium 
characterized by a coefficient 77 of viscosity and a shear modulus G, the so-called Maxwell model. In fact, 
the theory also takes into account effects of anisotropy, 47]. It predicts that in the region of hypersound, 
a few GHz or more, the attenuation coefficient for sound waves , a n , ceases to depend on frequency uj, 
whereas in the low frequency classical region, in which the Stokes law is valid, the dependence reads 
a n ~ uj 2 . Davis and VanZandt, (43j, used the approach of the Maxwell theory to find estimates for the 
damping of the DNA elastic modes, taking into account the first and the second hydration layers, and 
the quasi-crystallin structure of water in a neighborhood of the DNA. They found the attenuation to be 
two orders of magnitude smaller than that given by the Stokes law. 

In fact, there is a further reason for rejecting the approach based on the Stokes law. The water 
molecules form hydration shells of DNA, [31j, [32|]. The primary hydration shell comprises the water 
molecules immediately adjacent to the DNA, about 20 molecules per nucleotide pair, which constitute a 
medium different from bulk water. The secondary hydration shell is generally considered to be similar 
to bulk water. But, at the spatial scale of the diameter of the DNA-molecule, that is several tens A, 
the water in the second hydration shell, is hardly a condensed medium. Indeed, in this case one should 
have accommodated its local quasi-crystalline structure, described by the icosahedral model, [5^ ]. [53j . 
which is to result in sophisticated dynamical equations. The conclusion is that, presently, it is difficult, 
if possible at all, to construct accurate theoretical estimates for the attenuation of DNA-modes. 

From the experimental point of view the situation is more advanced. The DNA helical modes were 
observed in the experiments on the Raman. [l2| . [l3l ]. [3ll ] . and the far- infrared, (4||, scattering. Therefore, 
one may suggest that the attenuation effects due to viscosity should not preclude elastic modes of the 
DNA. At the same time small relaxation times for damping between the DNA and the first hydration 
layer, of order of several tens ps, (see [Uj]) should result in the double helix of DNA concerted motion 
with surrounding layer of water. The circumstance could be accommodated within the framework of the 
semi-phenomenological model of the present paper. In fact, the DNA molecule and its first hydration 
layer still form a helix structure, and the mutual motion of constituent bases of a pair together with 
hydration water molecules could be described with the field Y. Of course, the values of the model's 
constants, K, r, should be changed, and for the time being there is lack of information as to their size. 

It is also worth noting that the effects of dissipation in aqueous solutions, where a certain form of 
the Stokes law could be possible, and in films, or fibers, should be quite different. So far there has 
been no comprehensive theoretical analysis of the dissipation, which would allow for comparing the DNA 
dynamics in solutions and in films . Nonetheless, the interplay of internal vibration modes and sub- 
millimeter electromagnetic irradiation was registered in paper [461 ] , using Fourier transform spectroscopy 
and films of the double-stranded homopolymers poly[A]-poly[U] and poly[C]-poly[G] . Employing the 
concept of normal modes, or oscillators, of macromolecules, developed earlier for proteins, [541 ] . and used 
later for DNA, [55], Globus et al, [461 ] . made a numerical simulation of their experimental results, and thus 
obtained an estimate for the relaxational parameter 7, which has the meaning of oscillators dissipation. 
It turned out that in the range of frequencies several 10cm _1 , the best fit for 7 is less than 1cm -1 , 
depending on the conformation of an external electric field and a sample. This value of 7 is too large for 
Eq. ([3"5]) . but the region of frequencies studied in [46| is far from the edge of GHz region, so that one may 
consider the question of acceptable rate of dissipation as still open, and suggest that studying the effects 
of mw-radiation on the DNA modes may be instrumental for understanding the phenomenon. 
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We see that the elastic dynamics of the double helix could have enough structure for providing a 
means for stretching the hydrogen bonds of the base-pairs of DNA, or generating the HBS-modes. If 
the vibrational modes of the DNA are not overdamped by the ambient solvent, and the balance between 
energies supplied and dissipated is favourable, the maintenance of the HBS-modes could be expected 
at the edge of the HBS-zone. The best technique for studying the H-bond stretching still remains the 
Raman spectroscopy on which certain improvements have been made (see [56[ and references therein). 
Thus, the HBS-modes, and also the breathing modes, are well accessible from the experimental point of 
view. 

The choice of specific means for generating torsional excitations of the DNA is important and inter- 
esting. In this paper we have envisaged mw-irradiation of the DNA. In case the interpenetration of the 
acoustic and the HBS-modes takes place, mw-radiation could maintain the HBS-modes, if the power 
density is sufficiently large, 100 mW/cm 2 or more. It is important that there is no need for long ex- 
posures of the sample to the radiation. If the effect be sufficiently pronounced, it may result in the 
formation of the bubbles of broken H-bonds. At this point it is worth noting that our estimate for the 
critical power density, 100 mW/ cm 2 , is by orders of magnitude larger than that officially prescribed, i.e. 
0.2 - 0.1 mW/cm 2 . 



IV. PROTON TUNNELING INSIDE THE HYDROGEN BONDS OF THE DNA 



In the previous sections we have just considered a few specific situations which, nonetheless, indicate 
that the use of simple models and rough approximations is not sufficient for the study of conformations 
accessible for the DNA. The main conclusion is that a molecule of the DNA has an intrinsic structure 
that should be accounted for. Similar situations happen in continuum mechanics when it is necessary 
to consider a medium having an internal structure. But, the problem of the DNA is more sophisticated 
since the system under investigation is not a continuum medium, in fact, it is a macromolecule. Looking 
at the problem squarely, we have to acquiesce that we should work within the framework of nonlinear 
elasticity theory, if we wish to follow in this way. But the situation is still even more difficult owing to 
the necessity to allow for the intrinsic structure of 'the material'. Thus, at first sight the problem does 
not look tractable. 

It is reasonable to diminish the scope of regimes under the investigation and confine ourself mainly 
to problems that could be treated within the framework of the internal, or inter-strand, dynamics while 
assuming that deformations of the molecule of DNA are small on the mesoscale, that is a few persistence 
length. The assumption does not exclude the presence of external influences, as can be inferred from 
results of Section [TTJ There is another point we have to take into account: the dynamics of proton inside 
the hydrogen bonds between the base pairs. In contrast to the inter-strand modes of the DNA, they are 
essentially quantum modes. Therefore, we need a model that combines classical elasticity of the double 
helix and quantum dynamics of the protons. 

The model of this kind had been worked out by Davydov, I75| . for the needs of protein dynamics. It has 
been used for the dynamics of protons in the DNA in papers[8a| , [4l[ . It is important that the Davydov 
model has enough structure to accommodate both the inter-strand modes and the proton tunneling. The 
interaction between the two dynamics is a subject of great controversy, and whether it really takes place, 
or is merely speculative, depends on values of the elastic constants of the DNA, which are by no means 
precisely known. But if the numerical values are favourable, we may expect an interesting interplay 
between these forces. In what follows, we shall try to see what consequences could be inferred in this 
case. 

By changing both positions and mutual orientations of the base pairs, the inter-strand dynamics of the 
DNA should deform the hydrogen bonds between them, and thus have some bearing upon the protons. 
It should be noted that a proton effecting a hydrogen bond between two bases of the DNA, do not have 
a unique position of equilibrium. Under ordinary circumstances it occupies a position that corresponds 
to the bases being in the amino, or keto, forms for adenine- guanine, and cytosine-thymine, respectfully. 
The change in position of the proton result in the transition amino/ket — » imino/enol of the DNA base 
pairs. Does the transition influence elastic properties of the molecule of the DNA ? At this point we 
again come across the interplay between the microscopical and macroscopical dynamics of the DNA, and 
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FIG. 1: Pairing of Thymine- Adenine in the normal keto/amino forms in comparison with the tautomer-shifted 
enol and imino forms. 



its significance for the proton transport. The latter is of primary importance for the DNA, for among 
other things it could be a cause of spontaneous mutations. It also raises a question of whether irradiation 
with electromagnetic waves could result in generating the inter-strand modes, deforming the dynamics 
of proton tunneling, and causing genetic effects. 

Recall that according to the Watson-Crick hypothesis [la|. the double helix of the DNA molecule com- 
prises the two strands linked together by purine-pyrimidine base-pairs of adenine-thymine (AT) and 
guanine-cytosine (GC), the four chemicals A,T,G,C existing in various isomeric forms, or tautomers, that 
may change into one another (see FIG.[U for example). Under ordinary conditions the equilibrium shifts 
towards the ammo-form for adenine and guanine, and the keto-form for thymine and cytosine. But the 
imino-form for the adenine and cytosine, and the enol-form for guanine and thymine are also possible, 
even though rare; in fact, they correspond to concentrations of 10~ 4 to 10~ 5 moles/liter. The implica- 
tions wrought by the tautomeric transitions are important in that the sequence of base-pairs constitutes 
the genetic information of cell, so that exact copies of the DNA should be produced during the replication. 
In fact, the complimentarity between the bases may change if a tautomeric transition takes place, and 
other combinations become possible, 
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Genol * * T , G < > T eno i 

in contrast to the usual and stable ones 

A< >T , G< >C 

An opportunity for generating "unnatural" pairs arises also from the tunneling of protons in hydrogen 
bonds (see FIG. [T]), which results in the formation of the pairs 

\A < ► T) — r- (Ai m i no < > T eno i) (41) 

(G < ► C) s (G eno l < ► Cimino) 

During the replication, tautomeric transition driven by the proton tunneling in conjunction with the 
complimentarity according to (|4"Tj) may lead to the change of base-pairs 

(A< ►T) =► (G< >C) (42) 

(G < ► C) =► [A < ► T) 

and result in loss, or corruption, of genetic information, i.e. mutations, p^l [i~7| The specific case given 
by the diagram (|42|) is called transition mutations; it has the property of being reversible, i.e. able to go 
back to the wildlife type. 

The arguments given above constitute the main points of the theory of spontaneous mutations suggested 
by Crick and Watson. (l6l. Ii"7l. H8l 1571] It is based on the assumption that the transitory tautomeric shifts 
of base-pairs may occur during the replication, i.e. when two molecules of DNA are formed from a 
paired molecule, so that the double-stranded molecule is split into two single strands, each of which 
controls the synthesis of a new strand complimentary to itself with the help of the special enzyme called 
DNA polymerase. It has been realized that the latter plays an active role in the selection of bases 
at replication[58j], so that it may affect the mutation rates. Thus, tautomeric transitions are not a 
unique cause of mutation; the situation is more subtle, and many questions, of quite a classical nature, 
wait their solutions. Nonetheless, the original idea of Watson and Crick still conserves its appeal, and 
even more so as its new links with other phenomena related to the mutagenesis are brought to light 
(see[H,|60l). So, Robinson et aljg^, report that the enol tautomer of iG, that is 2'— deoxyisogine, may 
form at physiological temperature (37°) and pair with thymine in a Watson-Crick geometry; thus, iG 
being present as the nucleoside, results in the formation of incorrect base-pairs during in vitro replication. 
[6ll . H2, HH, [H, [66[ Robinson et al[6(|, suggests that iG ■ T pairing may have a bearing on mutagenesis 
in vivo involving tautomers of the common nucleobases. On the other hand, Fresco et al[65| . have found 
that the imino tautomer H0 5 dCyt may serve as an example of an unfavored base tautomer making for 
substituting mutagenesis. 

Mutations within the framework of the Crick- Watson model of DNA and in conjunction with the 
concept of tautomeric transition, have been drawing attention, beginning from the early fifties [lol [l6l 
E3, EE H3) to the present time, and involved the use of condensed matter theory. So, one of the first 
papers in this direction was published by Geracitano and Persico[67|, who suggested that there should be 
expected a collective behavior of codons, resembling that taking place in hydrogen-bonded ferroelectric 
crystals. 

In this paper we intend to look after the interplay between tautomeric transitions caused by the proton 
tunneling in base-pairs ( see FIG.[T]) and elastic properties of the double helix. 

We feel that quantum effects caused by the proton tunneling may have an appreciable bearing on 
mutagenesis. In this respect we would like to draw attention to the fact that mutations could be generated 
by irradiation with electromagnetic waves in infra-red region corresponding to the energies of tautomeric 
shifts in the base pairs of DNA. In fact, as was pointed out by Sukhorukov et al[68], the available data 
on the absorbtion spectra at v = 1697cm -1 , for the synthetic polynucleotide (PolyU)-(PolyA) forming a 
two-stranded structure, may indicate the transfer of protons between purine and pyrimidine bases. |71j In 
papers[68, 69] there is reported an absorbtion band at v — 1712cm -1 , for DNA at certain values of pH; 
the authors claim that the effect could be due to the proton tunnelling in hydrogen bonds of base pairs. 
Similar results are obtained in[7(| for two-stranded (PolyC) in LB-films. 

It is worth noting that the interplay among the proton tunneling and the elastic properties of DNA 
may manifest itself in the dynamical properties of mutagenesis. To be specific, it could result in the 
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FIG. 2: Mutation (black square) in a region different from a set of base pairs subject to an initial mutagenic 
action (arrows). 



phenomenon that the action imposed upon a set of base pairs of DNA may finally cause a substitution 
mutation in a different region of base pairs. One could expect the effect similar to the freak waves of 
nonlinear theory, when a low intensity initial perturbation for which the probability of mutation per base 
pair is low, could focus on a few base pairs and result in a mutation, ( see FIG. [5]). 

The main point about our hypothesis is that, since the 7r— electrons of the tautomeric rings of the 
nucleotides have direct bearing on the interaction of the plates of adjacent base-pairs[zl,[73j], a tautomeric 
transition of base-pairs should substantially influence the distribution of delocalized electrons of the 
nucleotides, i.e. the tt— electrons, and result in deformation of the elastic system of DNA. The hypothesis 
is in accordance with the conclusions of paper. [69( It is worth noting that tautomeric transitions may 
occur in several base pairs, not necessary adjacent, at a time, and their dynamics is determined by the 
proton tunneling. In fact, for one thing the latter is due to electrostatic interaction, i.e. the dipole 
forces, between the protons belonging to adjacent base-pairs, and for another the elastic system of the 
DNA molecule, which plays a role like that of the crystalline lattice of the polaron theory. According to 
D. Landau's original idea, a charge moving in the crystalline lattice, deforms the latter so that an effective 
field U e f f is generated, and within the framework of a self-consisted picture its motion is determined by 
U e ff. This argument can be also applied to the propagation of excitations of molecules constituting the 
lattice, that is the exciton theory. An important specific case is the regime in which the lattice relaxes 
to an equilibrium state fast enough so that we could neglect its motion and consider only the motion of 
the exciton, or charge. Then we may write down an effective hamiltonian that allows for the deformation 
of the lattice caused by the exciton, or charge, and next using an appropriate trial function derive an 
equation for its motion. [74| Davydov[75j], used the idea for the theory of a-helix in proteins, the molecule 
of protein playing the part of the crystalline lattice, and Volkenstein[73|, for his conformon theory. 

The essential point about the possible interplay among the proton tunneling and the conformation of 
the DNA is the values of its elastic constants. In fact, as far as the transport of torsional stress (torque) 
along DNA is concerned, its estimates obtained by various means diverge widely. The numerical values 
derived with the help of the theory of continuous media[8l|, are of the order t cx 10~ 17 dyne ■ cm [78l l80l|. 
whereas the experimental evidence (79j. indicates that it can attain the value of t a 10~ 13 dyne ■ cm. 
Philip Nelson [7g|, suggested that these deviations could be due to small bends in the helix backbone, so 
that one may assume 



For describing the elastic properties of the double- helix we may use the approach worked out in [82| , [83j , 
[85| . Thus, the double helix is considered as a one-dimensional lattice of vectors y n describing the mutual 
position of the two strands at sites corresponding to the base-pair of index n. It is important that the 
system has a twisted ground state characterized by the twist vector fi, so that the elastic energy of the 
molecule can be cast, at least for sufficiently small y n , in the form 



r cx 10 



,-17 



10 



-13 



dyne ■ cm 




(43) 



where the first term is the kinetic energy, the second one the elastic torsional energy and the last one 
corresponds to the separation of the two strands. The covariant derivative that accommodates the torsion 
of the molecule, reads 
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Here a is the spacing between the adjacent nucleotides, M is the mass of base-pair. For the sake of 
simplicity, we shall assume that the torsion vector fl is always parallel to the axis Oz, that is 

n = (0,0, n) 

and the vectors y n describe only transversal motions, that is y 3 = 0. It should be noted that we consider 
a very simplified model, use the harmonic approximation for its elastic energy, and assume that all 
sites, corresponding to base-pairs are identical. The subtle question is the value of the elastic constant 
K] obviously enough it has a direct bearing on the torque r mentioned above, and therefore, as was 
discussed above, its estimate may read 

K (x 10- 17 ^10~ 13 er ff 

It should be noted that the calculations within the framework of molecular dynamics, (see paper [26| and 
references therein), give the upper value for K, i.e. close to 10~ 12 4 10 -13 erg. 

The interplay between the torsional stress due to the relative motion of the base-pairs and the proton 
tunneling is very important. As was mentioned above the tautomeric transitions are driven by the proton 
tunneling, and therefore we shall describe them quantum mechanically, that is the stable amino/keto form 
corresponding to the ground state of proton, and the unstable imino/enol one to the excited state. [HI In 
accord with the qualitative character of our approach we neglect the fact that the tautomeric transitions 
in question, involve the tunneling of more than one proton, and assign only one proton to each site of 
the lattice. There are few hydrogen bonds in which the protons are transferred towards the imino/keto 
groups, or if one uses the concept of the two-level system, excited states. Therefore, one can consider the 
system as being close to equilibrium, or only weakly excited. This suggestion is very important for what 
follows. 

We shall describe the states of a base-pair at site n with the Bose operators b„ , b n that verify the usual 
conditions 

[ b n,Kn\ = $nm, [b n ,b m ] = [b+, 6+] = 0, 

and accommodate the assumption that the tunneling states of protons be described as two level systems 
by considering only their ground states and the first excited state. The energy of the protons, neglecting 
the interaction with the elastic degrees of freedom, reads [13] 

H P = Y,E b+b n + K £(6+& n+1 + b+ +1 6 n ) 

n n 

Here E is the energy of the tautomeric shift; its estimates depend on the choice of nucleotide and 
according to quantum chemistry calculations vary within the range of 2 4- 10 Kcal, (see[2j and references 
therein). The constant k could be ascribed to dipole interactions between adjacent sites, similarly to 
Davydov's theory. [75| Presently, there are no reliable estimates of its value (see below); by analogy 
with the Davydov theory one may assume that it should correspond to the characteristic frequency of 
tautomeric excitation due to the proton tunneling, that is of the order 10 11 Hz, or less. This estimate is 
generally accepted (see below). 

The central point of the model introduced in[84| is the interaction between the elastic degrees of freedom 
of DNA and the tautomeric transitions, or the proton tunneling in nucleotides; it reads 

H I = -\J2 (VfTn-£n) b+b n 

n 

Here the vectors h n give the spatial orientation for the hydrogen bonds of base-pairs 

h n = ( cos na, sinna, 0) 

It is important that the angle a is the rotation angle of the double- helix, and thus close to the angle Q, 
in accord with the fact that the covariant derivative term in the elastic energy provides only a qualitative 
description for the DNA-helix. 
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An argument in favor of the choice for the interaction Hi is that it takes into account the deformation 
of positions of adjacent base-pairs and its influence on the it— electrons of the bases, and therefore, 
the tautomeric transitions, or the related excitations of protons. According to the theory of [73| . the 
interaction could be appreciable. Thus, one may suggest that the interaction term could be larger than 
the tunneling term in the equation for Hp given above. 

Concluding we write the total energy within the framework of the model introduced in[84| in the form 



Htotal = H to r + Hp + Hj 

To find weakly excited states we shall use the Davydov approximation[75|,[86j], that is we shall look for 
the state vector of the system using the trial function 

\V>=J2 A n (t)-b+\0> 

n 

where |0 > is the ground state of the system, for which all the base-pairs, or the protons in the hydrogen 
bonds, are in the ground state, the amplitudes A n (t) being subject to the constraint 



The adiabatic approximation, which is important for the implementation of the Davydov theory, holds 
for the following reasons. The vectors y n describe the dynamics of base-pairs, that is relatively massive 
objects, and therefore one may consider them as classical fields. [8^. [8J|- We can derive the size of 
characteristic frequencies for y n from expression of the elastic energy. In fact, the mass M is that of 
the base-pair, that is of the order 500 Dalton, and K is of the same order of magnitude as the torque r 
discussed above. Hence, we get the characteristic velocity v for the y modes 




Interesting numerical values for the velocity v follow from the equation indicated above and the rough 
estimates for r or K we have mentioned. Indeed, for K cx 10~ 17 dyn ■ cm or less we obtain 

v oc 10 2 cm/sec 

For wavelengths of a few tens of A it gives the characteristic torsion or phonon frequencies of the order 

v y cx 10 8 10 9 Hz 

On the other hand, if we use the values for K provided by the molecular dynamics simulations [HJ , we get 
the velocity of excitations of the order 1000m/ sec, and v v cx 10 11 -j- 10 12 Hz, as for ordinary condensed 
media. 



The elasticity of the DNA strongly depends on nucleotide sequence, and therefore the arguments given 
above are only of qualitative nature. Coleman et al, [HI], put forward a lattice model of the DNA in which 
they try to accommodate the sequence dependence of elastic properties. At each side the deformation 
of the lattice is described by six kinematical variables: the three angular variables 9? (tilt, roll, twist) 
and the three displacement variables pf (shift, slide, rise). The elastic energy \& of a DNA segment is 
assumed to be the sum 

* = 5> B 

n 

of the interaction energies of adjacent base pairs tp n , which are functions of the above kinematic variables. 
In the notations of paper [HJ energy tp n reads 
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where Fg, H% are constants with Fgj = F£, H% = H£, i,j = 1,2,3. Coleman et al, [881, 

estimate 

F n = F 22 = 4.27 x 1(T 2 -f-s 

aeg z 

For T = 300 the above estimate gives 1.76 x 10~ erg, and the velocity of 'angular' waves corresponding 
to the above expression for the energy, several hundred raj sec. If the corresponding excitations of the 
double helix are of wavelength comparable with the distance between adjacent base pairs, we get the 
GHz-frcquency range. It should be noted that inter-strand modes, which are likely to correspond to the 
angular waves in the double helix, have been detected and measured in the region of tens - hundreds 
GHz, pjl, [H, Hi], H3 • Consequently, their propagation velocity could be well within a few hundred m/sec. 

It is instructive to compare the values of v y with the transition frequencies for tautomeric reactions 
inside the nucleotides, 



vp 



K 

2nh 



The estimates for the latter differ considerably, [871 ]. 

v P ot 10 6 10 u .Hz 

The lowest estimate, 10 6 Hz appears to be not unreasonable (V. Bendcrskii, and J.L.Leroy, personal 
communications) . 

The relative sizes of vp and u y are important for choosing the right approximation for the model. In 
fact, if we are at the lowest end of the spectra vp, then according to the estimate for v y obtained above 
the characteristic times for the acoustic modes are at least by an order of magnitude smaller than for the 
protons. In this case, we may suggest that the elastic system should follow the motion of the protons in 
hydrogen bonds, adjusting itself to it, so that a kind of adiabatic approximation can be employed. In 
this paper we shall follow this conjecture. 

Thus, we assume, as in paper [84j|. that the adiabatic approximation is valid, and therefore neglect the 
kinetic energy of the elastic system and take into account only its potential energy generated by the 
field y n . Then we are in a position to apply the self-consisted method of the exciton theory, in the form 
suggested by Davydov[75j], that is to calculate the mean value 

U eff =<V\H tor +Hj\V> (44) 

find the minimum, of U e ff with respect to y n , substitute it into the equation for the total energy 
Htotai so as to get the effective Davydov hamiltonian Tin, which depends only on the operator variables 
&+, b n , the classical variables y n having disappeared through the minimization. To make this paper as self 
consistent as possible we reproduce the necessary calculations of paper [84j]. Thus, we obtain an equation 
that has the form of the Schrodinger one 

ih^-\v>=n D \v> (45) 

at 

and in which the wave function \T> > should be of the form prescribed above. The assumption that the 
excited states correspond to the set of two-level systems is accommodated by the requirement that the 
operators are allowed only in the first power. It results in a system of equations, called the Davydov 
equations, for the amplitudes A n , which one obtains on equating the coefficients at 6+ on both sides of 
(|43| . ( see 75] for the details ). In this paper we consider the case of the stretching energy for the DNA 
strands being smaller than the torsional one, that is we assume 

ea 2 



being small enough. 
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The Davydov hamiltonian for our problem reads 

H D = Y< E ob+b n -J2<b+ +1 b n + b+b n+1 ) (46) 

n n 

n n 

A 2 ea 2 



^cosl m -"V • cos [(m - n)(0 - a)} \A m \ 2 \A n \ 
-, A - E cos |ro - n| • cos [(m - n)(0 - a)] |A n | 2 6+& m 

m.n 



ifO 2 
A 2 ea 



2 c „2 



and the equation for the amplitudes A„ 



ih-Q t A n = EaA n -K(A n+ i+A n -i)-—\A n \ 2 A n -—(J2\Am\ i )A n (47) 

m 

+ E c OS l mi -™^- COS [(m 1 -m 2 )(^-a)]|A mi | 2 |A m2 | 2 )A„ 



mi ,7712 

9 

ea 



if Kn 2 



(]T c OS l m -"l0 • cos [(m - n)(0 - a)] |A m | 2 )A„ 



which has the form of a nonlinear discrete Schrodinger equation for the amplitudes A n . The terms 
given by the third and the fourth lines of the equation written above, describe an interaction that has a 
very long range, so that one may claim Eq. (|47|) is essentially nonlocal in character. 

We shall make an important approximation that agrees with the general qualitative setting of our work, 
and set 



where 4> = arctan f2. Thus, the oscillating factors in Eq. (|47|) are cancelled out. With the help of the 
reduced variables B n 

A n = exp{- l -E t) B n (t) 
n 



we cast the equation for A„ in the form 



d A 2 
ih— B n = -n{B n+1 + B n -i) - —\B n \ 2 B n (48) 

^(£|£ ro | 4 )B n 



K 

in 

A 2 ea 2 
K KQ, 2 

7711 ,7712 



( cos lmi - m2[ <t>\B mi \ 2 \B m2 \ 2 )B n 



|^(E c - |ro - n| ^ s -i 2 )^ 



Introduce the characteristic frequencies 



k A K 

2nh 2-kTi 2irn 



(49) 
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and the dimensionless time 



T = t- vp 



It should be noted that the frequencies v y and vtor are not identical, v y ^ Vtor- Then the Davydov 
equation takes the form 

i^fBn = -{B n+1 +B n ^)-W\B n \ 2 B n (50) 
-W(J2\B m \ A )B n 



-WA( cos^- m ^4>\B mi \ 2 \B m2 \ 2 )B n 

mi ,m2 

-WA(J2cos^- n ^\B m \ 2 )B n 



in which 



W = (51) 

Vp ■ V t or 

A ■ <-» 

Now we aim at making the numerical simulation of Eq. (|50[) for various values of the parameters W, A, 
looking for solutions of the soliton type. We use the term soliton in a sense close to that used by applied 
scientists, i.e. a solution different from zero in a finite region of space, whose size we shall call the size of 
soliton, and preserving its shape for very long periods of time. For some values of W, A it has the form 
identical to the usual one, i.e. corresponding to the non-linear Schrodinger equation, but generally our 
solitons are different. The standard definition suggests that it be of the form 

Y(x,t) = e i( - qx -^ijj(x-vt) (53) 

in which ip is a real function. It is by no means clear that our solitons always have the form given by 
Eq.(53). 

The parameter A is a quantitative characteristic that enables us to take into account the structure of 
the double helix, and also the relative size of the torsional and deformation energies. In fact, A determines 
the magnitude of the nonlocal terms in Eq. (|50p , and in this respect it is worthwhile to note that for certain 
values of A and W we have not been able to find soliton solutions, e.g. A = 0.2 and W = 2, at least 
for physically reasonable sizes of solitons, i.e. less than 150 base pairs. The last constraint is due to the 
fact that we consider straight segments of DNA, parallel to Oz-axis, and therefore they should be of a 
size less than the persistence length, that is about 150 base pairs. But it is important that generally 
the condition A ^ does not forbid the existence of solitons, and its influence only results in the size of 
soliton becoming larger, which is quite natural, for A represents non-local terms in Eq. (|50[) . The general 
case of soliton with A not equal to zero, even though small, is illustrated in FIG. [3] To understand the 
general situation let us consider the two special cases. 

1. Stationary solutions in the sense that the absolute value, |B n (t)| does not depend on time. For 
the usual solitons given by Eq. (|50[) this requirement means that the velocity v — 0. The typical case 
is illustrated in FIG. Q] (a), for W — 10 and A = 0.5. The half-width of soliton equals to one spacing 
between base-pairs, that is the solution is extremely narrow, and according to our main hypothesis it 
must correspond to the tautomeric transition of a base-pair. The very interesting case is illustrated in 
FIG. [U (b), W = 5 and A = 0.5. There is a central peak of half- width 1.5 • a which stands still, and two 
symmetrical wave packets, moving in opposite outward directions. The distance traveled by these wave 
packets during 0.018 msec is equal to 33 base pairs. The value of vp was taken 10 6 Hz. 

2. The usual solitons given by Eq. (f53"|) . The half- width of these solitons may be several tens of base-pair 
spacings, and thus they could correspond to tautomeric transitions taking place in adjacent base-pairs. 
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FIG. 3: Typical moving solitons. Solid line: \A n \. Thin line: the real and the imaginary part of the amplitude 
A n , up — 10 6 Hz. (a) for W — 0.75, A = 0.001, velocity 1330 base pairs per msec, the period of time spent 
0.533 msec. (b) Moving soliton for W = 0.75, A = 0.075, velocity 1340 base pairs per msec, the period of time 
spent 0.530 msec, distance travelled 707 base pairs, (c) Moving soliton for W = 2, A = 0.1, velocity 850 base 
pairs per msec, the period of time spent 0.577 msec, distance travelled 491 base pairs. The values of W, A are 
close to the borderline [see FIG. [5] (b)], dividing the region of stable solitons from the unstable ones. 



The typical cases are illustrated in FIG.[3l It is interesting to note that these solitons move, even though 
slowly. Their velocity is given by the asymptotic formula 

v ~ 2 avp sin(ag) 

Hence, one might suggest the picture of tautomeric transitions moving along the DNA-moleculc. 
Both types of solutions indicated above are stable with respect to perturbation of W and A. 
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FIG. 4: Typical breathers. Solid line: \A n \. Thin lines indicate the real and the imaginary part of the amplitude 
An, vp — 10 6 Hz. (a) Breather, or still soliton. W = 10, A = 0.5. The period of time spent 0.501msec. (b) 
Radiation emitted from the motionless central peak during the period of time 0.018, for W = 5 and A = 0.5; the 
velocity of side waves 1830 base pairs/msec, distance traveled 33 base pairs. 



Perhaps, the most characteristic feature of discrete non-linear Schrodinger equation is solutions that 
periodically oscillate in time and decay exponentially in space, or breathers. From a purely qualitative 
point of view the existence of breathers can be inferred from a truncated version of Eq. (|50|) . Let us 
neglect all the terms on its RHS except the first two, that is consider 



th—^- = -(£?„+! + B n _x) - W\B n \ 2 B, 
at 

and look for B n such that 



Sivt 
n = e a, 

a n being real. Next, cast the equation for a n in the form 
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-(a„+i - 2a„ + a„_i) - [Wa 3 n + (2- e)]a„ = 

Suppose that the soliton we are looking for is large enough so that we may change the expression 
a n +i — 2a„ + a n -i for the second derivative. Thus we obtain the equation 



a" + [Wa 3 + (2 - e)]a = 
or the conservation law for one dimensional motion with the effective potential 



1/ 2 - £ 2 , W 4 

The soliton solution exists for e > 2, and its size tends to infinity as e — > 2. On the other hand for large 
W we may expect thin solitons. 

The key point is that the nonlocal terms generated by the double helix bring serious modifications to 
the picture given above. Thus, we may infer that the dimcnsionlcss constants W and A play a crucial 
role in determining the form of solitons for Eq. ([50|) . The general situation to the effect is illustrated in 
FIG. in which the horizontal axis corresponds to values of v, that is the soliton frequency measured in 
units of vp. It should be noted that v defines only the main Fourier component both for solitons and for 
breathers, so that v turns out to be only a rough characteristic. The breathers are represented by the 
solid line, that serves also as a right-handed border for the region of moving solitons. This line continues 
up to infinity with both W and v rising. It's lower end was not clearly found. The lower left-handed 
border of the soliton region is not strictly defined, owing to the fact that there are solitons for values of 
W and v lower than the borders but of sizes greater than 100 base pairs, that is beyond the physical 
context of our problem. The upper left-handed part of the border is determined by solitons turning out 
to be unstable for values of W and v beyond the boundary, and thus is not defined clearly. We see that 
the soliton region is decreasing as A grows (see FIG. [5]), and for A = 0.2, FIG. [5] (d), there are only 
breathers, at least under the constraint of their size being less than 100 base pairs. It is worth noting 
that Eg. ([50)1 derived in[84j is valid only for small A. 

We would like to draw attention to a class of solutions that correspond to the nomenclature of " freak 
waves" [8^, [9(3], and which may have a bearing upon the dynamics of tautomeric transitions. A solution of 
the type is illustrated in FIG. [SI It is characterized by an initial set of amplitudes B n (t) which is a broad 
distribution of the size of 80 base-pair spacings; after the period of time 0.017 msec, the characteristic 
frequency vp being taken 10 6 Hz, it focuses itself on a narrow peak of half-width of one spacing. The 
peak exists for the brief period of time 0.002 msec, and next breaks down into a broad distribution again, 
i.e. a kind of partial self focusing is taking place. Thus, there may exist low probability tautomeric 
transitions distributed over wide areas of the molecule, and which may collapse into a small region of 
the molecule, and stay there for a period of time, brief but perhaps sufficient to cause mutation. For 
finding the initial configurations producing the peaks indicated above we used the method of numerical 
integration backward in time, similar to that used in papers. [89 . |90| 

It is instructive to see the conformation of the field y n accompanying the dynamics of solitons. The 
typical configuration of y n corresponding to the soliton solution of Eg. ([50)) is shown in FIG.UJ 

We see that the formation and the dynamics of solitons corresponding to the tunneling of protons is 
accompanied by conformational changes. Thus, we may infer that the tautomeric transitions studied 
in this paper are to a large extent similar to conformons. [77| In fact, the concept of conformon was 
suggested for describing the dynamics of charge transfer in macromolecules, especially proteins, similarly 
to the situation we are considering. 

It is also interesting to invert the picture discussed above, and suggest that there is a deformation of 
the double helix described by a distribution of y n , like that shown in FIG. [71 then there should be a 
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FIG. 5: Sets of W and v that allow for soliton solutions. Solid line represents still solitons, or breathers. Shaded 
area represents moving solitons. Points are trial solutions. The transition from solitons to breathers goes smoothly 
as the soliton speed decreases to zero. We take into account only solutions of size less than 100 base pairs, (a) 
A = . (b) A = 0.1. (c) A = 0.15. (d) A = 0.2. 



distribution of the soliton amplitude, that is a proton tunneling generated by the conformational change. 
Thus, we may suggest that the conformational transitions may result in the proton tunneling, or the 
tautomeric shifts. 

As was shown above, the dynamics of tautomeric transitions in DNA depend on elastic properties of 
the latter and proton tunneling in base pairs; W and A serving as indicators for possible regimes. Our 
numerical simulation suggests that the interesting tautomeric dynamics may happen for W > 1. This 
allows for sufficiently wide range of material constants of DNA so that the phenomenon could occur. The 
second constant, A, provides a quantitative characteristic for the part played by the double helix; it can 
totally modify the structure of solitons corresponding to tautomeric transitions. 

Depending on the value of W one may expect the existence of two quite different dynamics: (1) solitons 
that move at velocities smaller by orders of magnitude compared with that of elastic excitations in DNA, 
and have a size of several tens of base-pairs, and (2) stationary solutions, or breathers, that have a form 
of peaks over a few base-pairs. We may suggest that the second type of solutions correspond to point 
mutations, whereas the first one may describe tautomeric transition moving along the chain of double 
helix, and therefore there may happen mutations related to the transition. Thus, one may suggest that an 
action imposed on a set of nucleotide in a region of the molecule might generate mutations in a different 
region owing to the motion of excitations corresponding to the proton tunneling. 

It is alleged to be known that by substituting the "artificial" nucleotides instead of the natural ones, 
e.g. brom-uracil for thymine, one can increase dramatically the rate of mutations; this could be due to the 
increase of tautomeric transitions inside base-pairs. At any rate, it is worthwhile to study the interplay 
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FIG. 6: Partial self-focusing of an initial low amplitude distribution on a peak for a period of time 0.002 msec, 
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FIG. 7: Moving soliton for W = 2 A = 0.1, velocity 880 base pairs per msec, up = 10 6 Hz . Maximum 
\y n \ = 0.52A (a) The amplitude A n . Solid line: \A n \. Thin line: the real and the imaginary part of the amplitude 
A n . (b) y„ distribution. Solid line: \y n \- Thin line: the first and second coordinates of the y n . 



between the rate of such transitions and mutations. Within the context of the present paper, artificial 
DNA of this kind could ease the stringent constraints imposed on W, as was indicated above. 

It is worth noting that the "focusing" of solutions (see FIG. [6]), similar to the freak waves which take 
place in the theory of non- linear waves [8^, [9fjj . may have a very important bearing on mutations. In 
fact, it amounts to the possibility of a weak external influence generating a low amplitude distribution of 
mutation sites that would focus itself later on a high amplitude distribution concentrated in a different 
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region of the molecule. Thus, one may expect generating mutations by low intensity agents distributed 
in a region of the molecule, or to put it the other way round, acting on a set of codons different from 
those that suffer the actual mutation. 
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